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The time-dependent Hartree-Fock calculation with a full Skyrme energy functional 
has been carried out on the three-dimensional Cartesian lattice space to study El 
giant dipole resonances (GDR) in light nuclei. The outgoing boundary condition 
for the continuum states is treated by the absorbing complex potential. The cal- 
culation for GDR in 16 O suggests a significant influence of the residual interaction 
associated with time-odd densities in the Skyrme functional. We also predict a 
O , large damping for superdeformed 14 Be at the neutron drip line. 



1. Time-dependent approach to nuclear response in the 
continuum 



d ' The quantum-mechanical problems are usually solved in the energy (time- 

independent) representation. Namely, we either solve an energy eigenvalue 
problem for bound states or, for scattering states, we calculate a wave 
function with a proper boundary condition at a given energy. However, 
if one wishes to calculate physical quantities in a wide energy region, the 
time-dependent approach is very useful because a single time propagation 
provides information for a certain range of energy. Another advantage may 
be an intuitive picture provided by the time evolution of the wave function. 

In Ref. 1, we have calculated molecular photoabsorption cross sec- 
tions in the electronic continuum by using the time-dependent and time- 
independent approaches. The results indicate the capability and efficiency 
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of the time-dependent method. In the present paper, we will show nu- 
clear response calculations using the same technique, the time-dependent 
Hartree-Fock (TDHF) method combined with the absorbing boundary con- 
dition (ABC). 

The TDHF state consists of A-occupied single-particle wave func- 
tions, {(f>i(r, i)}»=i,... t A, each of which is complex and has two components 
(spinor). The three-dimensional (3D) Cartesian coordinate, (x, y, z), is dis- 
cretized in rectangular lattice and derivatives are estimated with the nine- 
point formula. The time evolution is determined by 

&(r, t) = exp (-i J dt'h KF [p{t')^j <j>i{r, 0), (1) 

where the HF Hamiltonian, hnp[p(t)], depends on (pi(r,t) and their deriva- 
tives. The initial state, {<f>i(r, 0)}j = i i ... ,a, is chosen as the HF ground state 
wave function perturbed by an instantaneous external field, <pi(t — 0) — 
e lkF cf>®. Here, the coupling parameter, k, can be arbitrary but should be 
small to validate the Fourier analysis (linear response approximation) . The 
time variable is also discretized in a small step At and the exponential 
operator in Eq. (1) is approximated by the fourth-order expansion, 

exp ( -i dt'hn F [p(t')] \ » V -\ (h RF [p(At/2]) n . (2) 
V J " / n =o n - 

The time evolution of the expectation value of the external field, 
(^(t)\F\^(t)), is computed, then we utilize the Fourier transform to ob- 
tain the strength function, S(E) oc Ej2 n S(E - E n )\(n\F\0)\ 2 . 

Next, let us discuss the treatment of the continuum. During the time 
evolution of the TDHF state, a part of the wave function may escape from 
the nuclear binding and become an outgoing wave. In the time-dependent 
approach, it is very difficult to impose the outgoing boundary condition 
(OBC) explicitly, because the outgoing particles have different energies. 
Furthermore, the boundary condition for non-spherical systems requires 
involved computation. 1 Thus, instead of exactly treating the OBC, we use 
the absorbing boundary condition (ABC) which can be practically identical 
to the OBC. In the ABC method, the Cartesian coordinate space (r < 
i?o) is divided into two parts: an interacting region (r < R c ) and a non- 
interacting region (R c < r < Rq). We introduce a complex absorbing 
potential, —irj(r), active only in the non-interacting region. Then, we solve 
the TDHF equation of motion, Eq. (1), with the box boundary condition 
(BBC) at r = Rq, 0^(r, £)| r=jRo = for all occupied orbitals. This means 
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that the treatment of the continuum simply requires an addition of the 
absorbing potential to the HF Hamiltonian. 



/ihf [/>(*)] 



BBC at r=R 



(3) 



OBC at r=R c 

We have demonstrated accuracy and practicality of the ABC method both 
in the time-dependent 1 ' 2 ' 3,4 and the time-independent approaches. 5 ' 6 



2. Applications 

2.1. Giant dipole resonance in 16 O 

Giant dipole resonance (GDR) in 16 O has been studied extensively. 7 Since 
16 is a doubly-closed spherical nucleus, the continuum RPA 8 is applicable. 
A microscopic calculation with the Skyrme energy functional was done 9 
and showed a two-peak structure for GDR in 16 0. However, the calculation 
is not fully self-consistent because the residual Coulomb, spin-orbit, and 
a ■ a part of the interaction are neglected. According to the best of our 
knowledge, the continuum RPA calculation with the Skyrme functional has 
never been carried out fully self-consistently. Therefore, it is worth while to 
perform a fully self-consistent TDHF+ABC calculation and to investigate 
effects of the neglected part of the residual interaction. 

We use the Skyrme energy functional of Ref. 10 with the SGII parameter 
set. The 3D sphere of radius Rq = 22 fm is adopted as a model space. The 
mesh spacing is Ax w Ay « Az « 1 fm inside the interacting region 
but gradually increase up to A « 3 fm in the non-interacting region. After 
obtaining the HF ground state of 16 0, the isovector dipole field, F = z n — z p , 
is activated instantaneously at t — 0. The time evolution is computed up to 
T = 30 H/MeV. The result of the Fourier transform is shown in Fig. 1 (b). 
Apparently, we observe only one peak at E w 20 MeV, except for small 
peaks and a shoulder. This is inconsistent with results of the continuum 
RPA. The disagreement turns out to be due to the fact that the continuum 
RPA neglects part of the residual interaction. In the real-time calculation, 
we must use the same density functional as the one to define the HF ground 
state. Therefore, it is impossible to neglect the residual Coulomb and spin- 
orbit interaction. However, it is possible to neglect time-odd densities in the 
functional, because this does not affect the ground state. The calculation 
without the time-odd densities is shown in Fig. 1 (a). Now, the result 
becomes very similar to that of the continuum RPA and there are two 
peaks in the continuum. There are still some discrepancies related to the 
residual Coulomb and spin-orbit force. Our result indicates an important 
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role of the time-odd density part of the residual interaction. Since the 
time-odd part of the Skyrme energy functional is not well determined by 
the ground-state properties, this might provide a useful constraint. 
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Figure 1. Photoabsorption cross section for 16 O calculated with different Skyrme energy 
Junctionals, (a) Functional with SGII parameter set but neglecting all the time-odd 
densities, (b) Full Skyrme functional with the SGII parameter set. 



2.2. Giant dipole resonance in 8 > 14 £?e 

Utilizing the same parameters and model space as 16 O, we discuss El reso- 
nances in light deformed nuclei, 8 Be and 14 Be. The quadrupolc deformation 
of the HF ground states of these nuclei are as large as f3 ~ 1 (prolate) . Cal- 
culated El oscillator strengths are shown in Fig. 2. We can see a large 
deformation splitting of the GDR for 8 Be. The low-energy peak is sharper 
and located around 7 MeV of excitation energy. This peak is associated 
with the isovector dipole oscillation along the symmetry axis (z-axis). The 
higher peak is rather broad and located around 22 MeV. The difference 
in the broadening leads to a greater peak height for the low-energy res- 
onance, although the integrated oscillator strength is twice larger for the 
high-energy resonance. 

For 14 Be at the neutron drip line, the average peak positions for low- 
and high-energy resonances are similar to those in 8 Be. However, the excess 
neutrons significantly increase the peak width of both resonances. The large 
deformation splitting almost vanishes in the total strength (thick solid line) . 
It is worth noting that the present calculation takes account of the escape 
and the Landau damping width but does not include the spreading width. 
Some additional broadening might make the total strength look like a single 
peak with a large damping. 
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Figure 2. El Oscillator strength distribution for 8 Be and 14 Be calculated with the full 
Skyrme energy functional with the SGII parameters. The thin solid (dashed) line indi- 
cates the El oscillator strength associated with the oscillation parallel (perpendicular) 
to the symmetry axis. 

3. Conclusion 

We have carried out the Skyrme TDHF calculation with the ABC in real 
time for El resonances in the continuum. The time-odd density components 
in the Skyrme energy functional may influence considerably the isovector 
GDR strength distribution. The GDR in Be isotopes have been studied 
and the large deformation splitting and large broadening are predicted for 
a drip-line nucleus, 14 Be. 
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